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Summary 


Nomenclature 


A description  of  the  micro-air  vehicle  (MAY)  con- 
cept and  design  requirements  is  presented.  These  vehicles 
are  very  small  and  tlierefore  operate  at  chord  Reynolds 
numbers  below  200, 000  where  very  little  data  is  avail- 
able on  the  perfonnance  of  lifting  surfaces,  i.e.,  airfoils 
and  low  aspect-ratio  wings.  This  paper  presents  the  re- 
sults of  a continuing  study  of  tlie  metliods  tliat  can  be  used 
to  obtain  reliable  force  and  moment  data  on  thin  wings  in 
wind  and  water  tunnels.  To  this  end,  a new  platform  force 
and  moment  balance,  similar  to  an  already  existing  bal- 
ance, was  designed  and  built  to  perform  lift,  drag  and  mo- 
ment measurements  at  low  Reynolds  numbers.  Balance 
characteristics  and  validation  data  are  presented.  Results 
show  a good  agreement  between  published  data  and  data 
obtained  with  the  new  balance.  Results  for  lift,  drag  and 
pitching  moment  about  the  quarter  chord  with  the  exist- 
ing aerodynamic  balance  on  a series  of  thin  flat  plates 
and  cambered  plates  at  low  Reynolds  numbers  are  pre- 
smted.  They  show  diat  tlie  cambered  plates  offer  better 
aerodynamic  characteristics  and  performance.  Moreover, 
it  appears  tliat  the  trailing-edge  geometry  of  tlie  wings 
and  the  turbulence  intensity  up  to  about  1%  in  the  wind 
tumiel  do  not  have  a strong  effect  on  the  lift  and  drag 
for  diin  wings  at  low  Reynolds  numbers.  However,  the 
presence  of  two  endplates  for  two-dimensional  tests  and 
one  endplate  for  tlie  semi-infinite  tests  appears  to  have 
an  undesirable  influaice  on  tlie  lift  characteristics  at  low 
Reynolds  numbers.  The  drag  characteristics  for  thin  flat- 
plate  wings  of  aspect  ratio  greater  than  one  do  not  appear 
to  be  affected  by  tlie  endplates.  The  effect  of  tlie  endplates 
on  the  drag  characteristics  of  cambered-platc  wings  is  still 
under  investigation.  It  is  known,  however,  tliat  endplates 
do  have  an  effect  on  tlie  drag  and  lift  characteristics  of  a 
cambered  Eppler  61  airfoil/wiiig. 


Copyriglit  ©1999  by  Thomas  J.  Mueller.  Published  by  RTO  with 
pennission. 


Symbols 

AR 

Cd 

Cd 

Cl 

Cl 

Ci^  or  Ci^ 

Cf/Co 

^rnf4 

L/D 

M 

Rcc  or  Re 

Uoo 

a 

do 

b 

c 

sAR 

t 

a 

Ctstall 

r 

Subscripts 

max 

min 


full-span  aspect  ratio 

drag  coefficient  (3D) 

section  drag  coefficient  (2D) 

lift  coefficient  (3D) 

section  lift  coefficient  (2D) 

lift-curve  slope 

endurance  parameter 

pitcliing  moment  coefficient  about 

the  quarter  chord 

slope  of  pitching  moment  curve 

lift-to-drag  ratio 

resolution  of  A/D  converter 

root-chord  Reynolds  number 

freestream  velocity 

lift-curve  slope 

2D  lift-curve  slope 

wing  span 

root-chord  length 

quantization  error 

semi-span  aspect  ratio 

wing  tliickness 

angle  of  attack 

zero-lift  angle  of  attack 

stall  angle  of  attack 

Glauert  parameter 


maximum 

minimum 


Abbreviations 


2D 

3D 

A/D 

TE 

UND-FBl 

UND-FB2 


two-dimensional  (airfoil) 
tliree-dimensional  (wing) 
analog-to-digital 
trailing  edge 

old  Notre  Dame  aerodynamic  force  balance 
new  Notre  Dame  aerodynamic  force  balance 
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Introduction 

There  is  a serious  effort  to  design  aircraft  that  are 
as  small  as  possible  for  special,  limited-duration  military 
and  civil  missions.  These  aircraft,  called  micro-air  ve- 
hicles (MAVs)  (Davis  et  al,  1996;  Ashley,  1998;  Wilson 
1998;  Dornheim,  1998;  Mraz,  1998;  andFulghum,  1998), 
are  of  interest  because  electronic  surveillance  and  detec- 
tion sensor  equipment  can  no^v  be  miniaturized  so  that 
the  entire  payload  mass  is  about  18  grams.  The  advan- 
tages of  a MAV  include  compact  system  transportable  by 
a single  operator,  rapid  deployment,  real-time  data,  low 
radar  cross-section,  difficult  to  see  and  very  quiet.  The 
potential  for  low  production  cost  is  also  an  advantage. 
The  primary  missions  of  interest  for  fixed  wing  MAVs 
include  surveillance,  detection,  communications,  and  the 
placement  of  unattended  sensors.  Surveillance  missions 
include  video  (day  and  night)  and  infrared  images  of  bat- 
tlefields (referred  to  as  tlie  ‘river  the  hill”  problem)  and 
urban  areas  (referred  to  as  “around  the  comer”).  These 
real-time  images  can  give  the  number  and  location  of  op- 
posing forces.  This  type  of  information  can  also  be  useful 
in  hostage  rescue  and  countar-dmg  operations.  Because  of 
the  availability  of  very  small  sensors,  detection  missions 
include  the  sensing  of  biological  agents,  chemical  com- 
pounds and  nuclear  materials  (i.e.,  radioactivity).  MAVs 
may  also  be  used  to  improve  communications  in  urban 
or  otlier  environments  where  full-time  line  of  sight  oper- 
ations are  important.  The  placement  of  acoustic  sensors 
on  the  outside  of  a building  during  a hostage  rescue  or 
counter-drug  operation  is  another  possible  mission. 

The  requirements  for  fixed  wing  MAV s cover  a wide 
range  of  possible  operational  environments  including  ur- 
ban, jungle,  desert,  maritime,  mountains  and  arctic  envi- 
ronments. Furthermore,  MAVs  must  be  able  to  perform 
their  missions  in  all  weather  conditions  (i.e.,  precipita- 
tion, wind  shear,  and  gusts).  Because  these  vehicles  fly 
at  relatively  low  altitudes  (i.e.,  less  than  100  m)  where 
buildings,  trees,  hills,  etc.  may  be  present,  a collision 
avoidance  system  is  also  required. 

The  long  term  goal  of  this  project  is  to  develop  air- 
craft systems  with  a mass  of  less  than  30  grams,  have 
about  an  eight  centimeter  wing  span  that  can  fly  for  20 
to  30  minutes  at  between  30  and  65  km/hr.  The  cur- 
rent goal  is  to  develop  aircraft  with  a 15  centimeter  wing 
span  that  have  a mass  of  about  90  grams.  The  gross 
mass  of  micro-air  vehicles  and  other  flying  objects  ver- 
sus Reynolds  number  is  shown  in  Figure  1,  with  the  data 
from  Jackson  (1996-97),  Taylor  (1969-70),  and  Tennekes 
(1996).  Since  it  is  not  possible  to  meet  all  of  the  design 
requiremmts  for  a micro-air  vehicle  with  currmt  technol- 
ogy, research  is  proceeding  on  all  of  the  system  compo- 
nents at  various  government  laboratories,  companies  and 
universities. 

Design  aims 

The  design  requirements  cover  a wide  range  when 
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Figure  I;  Reynolds  number  range  for  flight  vehicles 


one  considers  the  diversity  of  possible  applications  for 
micro-air  vehicles.  The  MAV  must  be  designed  as  a 
system  consisting  of  airframe,  propulsion,  payload  and 
avionics.  Although  much  smaller  than  currently  opera- 
tional UAVs,  electrically  powered  MAVs  will  have  ap- 
proximately the  same  weight  fractions,  that  is,  21%  for 
the  airframe,  11%  for  the  engine,  30%  for  the  battery,  21% 
for  the  payload,  and  17%  for  avionics  and  miscellaneous 
items.  Minimum  wing  area  for  ease  of  packaging  and 
pre-launch  handling  is  also  important.  Figure  2 presents 
the  payload  mass  versus  wingspan  for  MAVs  and  otlier 
larger  UAVs. 
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Figure  2:  UAV payload  wingspan  (Davis,  1999) 

(Reprinted  with  permission  of  MIT  Lincoln 
Laboratory,  Lexington,  Massachusetts) 


A typical  fixed  wing  MAV  mission  (Morris,  1997) 
could  include  the  following  sequence  of  events: 

1.  Launch  and  climb  to  100  meters 

2.  High  speed  dash  (64 /cm/hr  Indicated  Air  Speed)  to 
target  (at  40  km  (hr  head  wind) 
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3.  Loiter  over  target  area 

4.  Maneuver  over  target  during  loiter  while  turning  at 
the  minimimi  radius 

5.  Descend  and  climb  over  target  area 

6.  Climb  to  100  meters 

7.  High  speed  dash  (64 /c/ny/rr  Indicated  Air  Speed)  to 
launch  point  (tail  wind  40  km/ hr). 

Mission  constraints  in  this  simulation  include  dura- 
tion, operational  radius,  minimum  turning  radius,  min- 
imum climb  angle,  maximum  altitude  and  number  of 
climbs.  Several  MAV  designs  have  been  built  and  flown 
witli  tliis  type  of  mission  in  mind.  A 15c?77  , square  plan- 
fonn  internal  combustion  engine  powered  vehicle  called 
the  Flyswatter  has  been  flown  by  Morris  (1997).  A nidder 
and  elevator  surfaces  arc  used  to  control  this  MAV.  The 
first  electric  powered  15  cm  MAV  witli  proportional  radio 
control  carrying  a video  camera  was  designed  and  flown 
by  Mattlicw  T.  Keainon  of  AewVnxmment.  This  vehicle 
called  the  Black  Widow  currently  holds  tlie  record  for 
endurance  at  22  minutes  (Keennon,  1999).  Other  MAVs 
witli  larger  dimensions  have  been  designed  and  flown  to 
help  develop  the  electronic  packages  and  control  systems 
(Harris,  1999;  and  Ailinger,  1999).  Although  these  are 
examples  of  current  vehicles,  further  improvements  will 
be  made  when  more  data  on  low  Reynolds  number  aero- 
dynamics is  available  and  smaller,  more  efficient  electric 
motors  and  propellers  have  been  developed. 

The  airfoil  section  and  wing  planform  of  tlie  lifting 
surface  occupy  a central  position  in  all  design  procedures 
for  flying  vehicles.  Therefore,  all  low  Reyitolds  num- 
ber vehicles  share  tlie  ultimate  goal  of  a stable  and  con- 
trollable vehicle  with  maximum  aerodynamic  cfTiciency. 
Aerodynamic  cfTiciency  is  defined  in  terms  of  the  lifl-to- 
drag  ratio.  Airfoil  section  and  {Ci/Cd)^^^^^^ 

as  a function  of  Reynolds  number  arc  shown  in  Fig- 
ures 3a,  3b,  and  3c  after  McMasters  and  Henderson 
(1980).  It  is  clear  from  this  figure  tliat  airfoil  perfor- 
mance deteriorates  rapidly  as  the  chord  Reynolds  num- 
ber decreases  below  100, 000.  Wliile  tlie  maximum  lift- 
to-drag  ratio  for  most  low-speed  fixed-wing  aircraft  ( 
T^oo  < 50  m/s)  is  greater  tlian  10,  values  for  insects 
and  small  birds  arc  usually  less  than  10.  Furthermore,  to 
achieve  tliese  values  for  MAVs  at  low  Reynolds  numbers, 
the  wings  must  emulate  bird  and  insect  wings  and  be  very 
thin  (i.e.,  i/c  < 0.06)  with  a modest  amount  of  camber. 

Requirements  for  a typical  propeller  driven  MAV, 
for  example,  include  long  fliglit  duration  (i.e.,  high  value 
of  C}y^/Cj9  at  speeds  up  to  65krn/hr  at  chord  Reynolds 
numbers  from  about  45,000  to  180,000  and  altitudes 
from  30  to  100  meters).  Since  these  vehicles  are  es- 
sentially small  flying  wings,  there  is  a need  to  develop 
efficient  low  Reynolds  number,  low  aspect-ratio  wings 
which  arc  not  overly  sensitive  to  wind  shear,  gusts,  and 
the  roughness  produced  by  precipitation.  Furthemiore, 


REYNOLDS  NUMBER  ~Rn 


(a)  Maximum  lift  coefficient 


REYNOLDS  NUMBER  ^ Rn 


(b)  Minimum  drag  coefficient 


(c)  Maximum  lift-to-drag  ratio 

Figure  3:  Airfoil  performance 
(McMasters  and  Hendersori  1980) 
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confidence  that  the  operational  vehicle  will  perform  as 
designed  is  important  in  all  applications. 

Flow  problems 

Although  design  methods  developed  over  the  past 
35  years  produce  efficient  airfoils  for  chord  Reynolds 
numbers  greater  than  about  200, 000,  these  methods  are 
generally  inadequate  for  chord  Reynolds  numbers  below 

200, 000,  especially  for  very  thin  airfoils.  In  relation  to 
the  airfoil  boundary  layer,  important  areas  of  concern  are 
the  separated  regions  which  occur  near  the  leading  and/or 
trailing  edges  and  transition  from  laminar  to  turbulent 
flow  if  it  occurs.  It  is  well  known  that  separation  and 
transition  are  highly  sensitive  to  Reynolds  number,  pres- 
sure gradient,  and  the  disturbance  environment.  Transi- 
tion and  separation  play  a critical  role  in  determining  the 
development  of  the  boundary  layer  which,  in  turn,  affects 
the  overall  performance  of  the  airfoil.  The  aerodynamic 
characteristics  of  the  wing  and  other  components  in  turn 
affect  the  static,  dynamic  and  aeroelastic  stability  of  the 
entire  vehicle.  Therefore  the  successful  management  of 
the  sensitive  boundary  layer  for  a particular  low  Reynolds 
number  vehicle  design  is  critical. 

The  survey  of  low  Reynolds  numb  O'  airfoils  by 
Carmichael  (1981),  althougli  almost  two  decades  old,  is 
a very  useful  starting  point  in  the  description  of  the  char- 
acter of  the  flow  over  airfoils  over  the  range  of  Reynolds 
numbers  of  interest  here.  The  following  discussion  of 
flow  regimes  from  1,000  < Rcc  < 200,000  is  a modi- 
fied version  of  CarmichaeTs  original  work. 

• In  the  range  between  1,000  < Re^  < 10,000,  the 
boundary  layer  flow  is  laminar  and  it  is  very  difficult 
to  cause  transition  to  turbulent  flow.  The  dragon  fly 
and  the  house  fly  are  among  the  insects  that  fly  in 
this  regime.  The  dragon  fly  wing  has  a sawtooth 
single  surface  airfoil.  It  has  been  speculated  that 
eddies  in  the  trouts  help  keep  the  flow  from  sepa- 
rating. The  house  fly  wing  has  large  numbers  of  fine 
hair-like  elements  projecting  normal  to  the  surface. 
It  is  speculated  that  these  promote  eddy-induced  en- 
ergy transfer  to  prevent  separation.  Indoor  rubber- 
powered  type  model  airplanes  also  fly  in  this  regime. 
It  has  been  found  that  both  blunt  leading  and  trailing 
edges  enhance  the  aerodynamic  performance. 

• For  chord  Reynolds  numbers  between  10, 000  and 

30,000,  tlie  boundary  layer  is  completely  laminar 
and  artificial  tripping  has  not  been  successful.  Ex- 
perience with  hand-launched  glider  models  indicates 
that  when  the  boundary  layer  separates  it  does  not 
reattach. 

• The  range  30, 000  < Re^  < 70, 000  is  of  great  in- 
terest to  MAV  designers  as  well  as  model  aircraft 
builders.  The  choice  of  an  airfoil  section  is  very  im- 
portant in  this  regime  since  relatively  thick  airfoils 
(i.e.,  6%  and  above)  can  have  significant  hysteresis 


effects  caused  by  laminar  separation  with  transition 
to  turbulent  flow.  Also  below  chord  Reynolds  num- 
bers of  about  50, 000,  the  free  shear  layer  after  lami- 
nar separation  normally  does  not  transition  to  turbu- 
lent flow  in  time  to  reattach.  Near  the  upper  end  of 
this  range,  the  critical  Reynolds  number  can  be  de- 
creased by  using  boundary  layer  trips.  Thin  airfoil 
sections  (i.e.,  less  than  6%  thick)  at  the  upper  end 
of  this  regime  can  exhibit  reasonable  performance. 

• At  Reynolds  numbers  above  70,000  and  below 

200, 000,  extensive  laminar  flow  can  be  obtained  and 
therefore  airfoil  performance  improves  although  the 
laminar  separation  bubble  may  still  present  a prob- 
lem for  a particular  airfoil.  Small  radio  controlled 
model  airplanes  fly  in  this  range. 

• Above  Rcc  of  200,000,  airfoil  performance  im- 
proves significantly  and  there  is  a great  deal  of  ex- 
perience available  from  large  soaring  birds,  large  ra- 
dio controlled  model  airplanes,  human  powered  air- 
planes, etc. 

Laminar  separation  bubbles  occur  on  the  upper  sm- 
face  of  most  airfoils  at  Reynolds  numbers  above  about 

50, 000.  These  bubbles  become  larger  as  the  Reynolds 
number  decreases,  usually  resulting  in  a rapid  deteriora- 
tion in  performance,  i.e.,  substantial  decrease  in  L/D. 
In  principle  the  laminar  separation  bubble  and  transition 
can  be  artificially  controlled  by  adding  the  proper  type  of 
disturbance  at  the  proper  location  on  the  airfoil.  Wires, 
tape  strips,  grooves,  steps,  grit,  or  bleed-through  holes  in 
the  airfoil  surface  have  all  been  used  to  have  a positive 
influence  on  the  boundary  layer  in  this  critical  Reynolds 
number  region.  The  type  and  location  of  these  so-called 
“turbulators”  and  their  actual  effect  on  the  airfoil  bound- 
ary layer  has  not  been  well  documented.  Furthermore,  the 
addition  of  a turbulator  does  not  always  improve  the  air- 
foil performance.  In  fact,  how  the  disturbances  produced 
by  a given  type  of  turbulator  influence  transition  is  not 
completely  understood. 

As  a result  of  this  critical  boundary  layer  behavior, 
several  important  questions  must  be  addressed: 

1.  What  is  the  free  stream  disturbance  level  and  flight 
envirorunent  for  a given  low  Reynolds  number  ap- 
plication? 

2.  If  the  flight  conditions  are  known  and  a suitable 
design  technique  was  available,  could  the  resulting 
vehicle  or  component  be  adequately  evaluated  in  a 
wind  turmel  which,  in  general,  has  a different  dis- 
turbance level  and  environment  than  the  flight  con- 
dition? 

3.  Is  the  hysteresis  in  aerodynamic  forces  observed  in 
low  turbulence  wind  turmel  experiments  present  in 
powered  applications  (i.e.,  do  structural  vibrations 
originating  with  the  propulsion  or  drive  system  affect 
boundary  layer  transition)? 
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4.  Because  the  critical  quantities  measured  in  wind  tun- 
nel experiments  are  very  small,  what  is  the  level 
of  accuracy  needed  to  improve  design  and  analysis 
methods? 

Preliminary  experiments 

Many  of  the  problems  plaguing  very  low  Reynolds 
number  research  involve  tlic  difficulties  associated  with 
making  accurate  wind/water  tunnel  models  and  obtain- 
ing reliable  data.  Because  the  boundary  layers  are  sen- 
sitive to  small  disturbances,  accurate  wind/watcr  tunnel 
models  are  very  important  in  the  evaluation  of  a given 
design.  Furthcmiorc,  because  the  forces,  pressure  differ- 
ences and  velocities  are  extremely  small,  a great  deal  of 
care  must  be  exercised  to  obtain  accurate  and  meaningful 
data.  Low  Reynolds  number  aerodynamic  research  has 
been  in  progress  at  tlie  University  of  Notre  Dame  since 
1978.  However,  chord  Reynolds  numbers  below  about 
80, 000  were  seldom  of  interest  in  tlie  studies  before  1996. 
Also,  most  of  the  studies  were  for  relatively  thick  airfoils, 
e.g.  the  11%  thick  Lissaman  7769,  the  13%  tliick  Miley 
M06-13-128,  and  tlie  13%  tliick  Wortmaim  FX  63-137 
airfoils.  The  only  relatively  thin  airfoils  studied  were  the 
Epplcr  61  and  Pfenninger  048  airfoils.  The  Eppier  61, 
shown  in  Figure  4,  was  originally  designed  for  model  air- 
planes witli  a chord  Reynolds  numbers  of  about  80, 000 
and  has  a tliickness  of  5.63%  and  6.3%  camber.  Figure  5 
shows  a schematic  of  the  Pfciminger  048  airfoil  geome- 
try tested  by  Bums  (see  Bums,  1981).  This  airfoil  has  a 
thickness-to-chord  ratio  of  4.8%  and  a 4.2%  camber. 


Figure  4:  Eppier  61  airfoil  profile 


Figure  5:  Pfenninger  048  airfoil  profile 


The  wind  tumiel  data  shown  in  Figure  6 for  the 
Eppier  61  and  the  Pfenninger  048  airfoils  was  obtained 
in  1980  and  published  by  Bums  (1981)  and  Mueller  and 
Bums  (1982).  Figure  6 indicates  that  for  chord  Reynolds 
numbers  below  90, 000,  tlie  tliinner  and  sharper  leading 
edge  Pfenninger  airfoil  perfonns  better  than  tlie  Eppier  61 
airfoil. 


Figure  6:  Maximum  lift' to- drag  ratio  versus  chord 

Reynolds  number  for  the  two-dimensional  Eppier  61  and 
Pfenninger  048  airfoils  (Burns,  1981) 


Bums  also  studied  the  flowfield  over  the  Eppier  61 
airfoil  for  differoit  Reynolds  numbers  using  the  smoke 
wire  teclmique.  The  location  of  tlie  boundary  layer  sepa- 
ration could  be  obtained  from  his  flow  visualization  pho- 
tographs. Figure  7 shows  tlie  effect  of  changing  the  angle 
of  attack  on  Uie  boundaiy  layer  for  Rcc  ==  46, 000. 

Flow  visualization  was  also  conducted  on  the  Pfen- 
ninger 048  airfoil  at  different  Reynolds  numbers  and  dif- 
ferent angles  of  attack.  Figure  8 shows  examples  of  the 
flow  visualization  images  for  Re^  ^ 47, 000. 

A new  series  of  experiments  was  perfomied  in  the 
Spring  of  1997  to  evaluate  several  thin  airfoil  shapes  us- 
ing the  existing  strain  gauge  force  balance  (UND-FBl)  in 
the  Hessert  Center  water  tunnel.  The  results  for  lift  and 
drag  from  these  experiments  down  to  a chord  Reynolds 
of  25, 000  were  very  encomaging.  A more  complete  ex- 
perimental study  of  one  of  these  airfoil  shapes  (i.e.,  the 
Eppier  61)  was  perfonned  during  the  summer  of  1997 
(Prazak  and  Mueller,  1997).  These  experiments  covered 
the  Reynolds  number  range  from  12, 000  to  63, 000.  Hy- 
drogen bubble  flow  visualization  was  used  to  detennine 
the  location  of  boimdary  layer  separation  and  the  existing 
wind  tunnel  lift/drag  force  balance  (UND-FBl)  was  used 
to  make  aerodynamic  measurements.  All  of  tliese  2D  ex- 
periments included  two  endplates.  The  (Q/C^i)max  for 
the  2D  Eppier  61  from  these  water  tunnel  experiments  is 
included  in  Figure  6 (A).  The  results  of  tlie  Prazak  and 
Mueller  study  in  the  water  tunnel  indicated  diat  by  adding 
a computer  data  acquisition  system  to  the  UND-FBl  force 
balance,  tlie  uncertainty  in  tlie  force  measurements  down 
to  Rcc  = 20,000  could  be  reduced  significantly  for  the 
full  span  models. 

Scope  of  present  study 

The  purpose  of  the  present  work  is  to  present  and 
discuss  tlie  measurement  problems  associated  with  small 
aspect  ratio  wings  at  Reynolds  numbers  below  200, 000. 
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Figure  8:  Smoke  flow  visualization  on  Pfenninger  048 
Figure  7:  Smoke  flow  visualization  on  Eppler  61  airfoil  airfoil  at  Rcc  — 47,000  (Burns,  1981) 

at  Rcc  — 46,  000  (Burns,  1981) 
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Both  wind  tunnel  cUid  water  tunnel  experiments  were  per- 
fonned  in  mi  attempt  to  acquire  2D  airfoil  mid  finite  low 
aspect  ratio  wing  data.  vSludies  of  tlie  effect  of  two  end- 
plates  on  Uie  results  for  2D  configurations  and  one  end- 
plate  for  semi-spmi  or  half  models  were  iilso  made. 

Apparatus  and  procedures 

Wind  tunnel 

The  Hessert  Center  for  Aerospace  Research  is 
equipped  witli  two  similar,  horizontal,  subsonic  open- 
circuit  wind  tunnels.  Each  indr^ift  tunnel  has  a contrac- 
tion ratio  of  20.6:1.  The  cross-sections  of  tlie  enlrmice 
and  test  section  tire  square.  The  hirgest  test  section  is 
two  feet  by  two  feet  (61  cm  by  61  cm).  The  contraction 
cones  are  designed  to  provide  very  low  turbulence  levels 
in  tlie  test  section.  Just  tiliead  of  tlie  contraction  cone 
tire  twelve  anti- turbulence  screens.  Botli  tlie  contraction 
cone  mid  tlie  test  sections  tire  mounted  on  rollers  to  pro- 
vide tin  easy  mcmis  of  interchanging  tliese  components. 
DownsLremn  of  tlie  test  section  is  tlie  diffuser  which  is 
fixed  into  the  wall  of  Uie  laboratory.  The  diffuser  de- 
celerates tlie  tiir  mid  also  gradually  transfomis  tlie  squtire 
contour  to  a circle.  The  impeller  is  driven  by  a variable 
speed  electric  motor.  By  vtirying  Uie  speed  of  Uie  motor, 
Uie  tunnel  speed  may  reach  a maximum  of  approximately 
12G/f/6‘cc  (36.6m/s)  wiUi  a four  square  foot  test  section 
(3, 221  cm^).  The  test  sections  used  are  six  feet  (1.82  me- 
ters) in  length.  All  wind  tunnel  experiments  presented  in 
this  report  were  conducted  in  one  of  Uiese  subsonic  wind 
tunnels.  The  range  of  velocities  required  for  tests  up  to 
Rcc  — 200, 000  could  easily  be  obtained.  In  general,  Uie 
minimum  velocity  for  force  balance  measurements  was 
kept  above  5 m/.s  (16.4  / f/s).  Figure  9 is  a schematic  of 
the  wind  tunnel  used.  The  freestream  turbulence  intensity 
was  approximately  0.05%  over  Uie  range  of  interest. 


Figure  9:  Schematic  of  the  low-speed  wind  tunnel 
(all  dimensions  in  mm) 


Water  tunnel 

For  tests  between  20,000  < Rcc  < 80,000,  an 
Eidetics®  free-surface  water  tunnel  wiUi  a 15  m x 18  in 


(38.1  cm  X 45.7  cm)  test  section,  pictured  in  Figure  10, 
was  used.  The  water  tunnel  is  also  located  in  Uie  mmn 
laboratory  of  Uie  Hessert  Center.  Water  velocities  up  to 
1.28  ft/ sec  (39  cm/ sec  ) can  be  obtained  in  tlie  test 
section.  A freestream  turbulence  intensity  of  less  Uiaii 
1%  has  been  reported  by  Uie  manutheturer  of  the  water 
tunnel. 


7.  DISCHARGE  PLENUTvl 

8.  RETORN  PIPING  (I^IEVERSIBLE) 

9.  FILTER  SYSTEM 

Figure  10:  Schematic  of  the  flow  visualization  water 
tunnel  (Eidetics®) 


Flow  visualization  techniques 

A qualitaUve  visucil  examination  of  Uie  flow  as  it 
passes  mi  airfoil  is  key  in  understmiding  its  qumititative 
aerodyncimic  characteristics.  In  Uie  wind  tunnel,  two  dif- 
ferent methods  can  be  used  to  generate  smoke  for  this 
flow  visualization.  WiUi  Uie  first  method,  smoke  is  gen- 
erated by  a device  which  allows  kerosene  to  drip  onto 
electrically  heated  filaments;  Uie  smoke  is  Uien  funneled 
to  a smoke  rake.  The  rake  has  a filter  bag  and  cool- 
ing coils  which  reduce  Uie  smoke  temperature  to  approxi- 
mately ambient  before  passing  Uirough  Uie  anti-turbulence 
screens  and  into  Uie  test  section.  With  Uie  second  meUiod, 
a fine  wire  is  placed  upstiemn  of  Uie  model.  This  wire 
is  coated  wiUi  oil  and  an  electric  current  is  applied  to 
Uie  wire.  As  Uie  wire  gets  heated,  Uie  small  beads  of  oil 
fonned  on  Uie  wire  burn,  which  gives  rise  to  fine  smoke 
slreaklines.  This  technique,  which  has  been  described  in 
more  details  by  Batill  mid  Mueller  (1980),  is  rel^rred  to  as 
the  smoke-wire  technique.  The  water  tunnel  is  excellent 
for  flow  visualization  using  eiUier  Uie  hydrogen  bubble 
or  dye  injection  technique.  A Kodak  DC  120  digital  emn- 
era  mid  CCD  video  cameras  are  available  to  capture  flow 
visualization  results. 
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Description  of  the  UND-FB1  balance 

Most  of  tlie  results  on  thin  plates  were  obtained  with 
the  existing  three-component  platform  aerodynamic  bal- 
ance UND-FB 1 . This  balance  can  be  used  to  measure  lift, 
drag  and  pitching  moment  about  the  vertical  axis.  The 
balance  is  an  external  balance  placed  on  top  of  the  test 
section  of  either  of  the  two  low-speed  wind  tunnels.  With 
this  balance,  lift  and  drag  forces  are  transmitted  through 
the  sting  which  is  moimted  directly  to  the  moment  sensor. 
The  moment  sensor  is  rigidly  mounted  to  the  adjustable 
angle  of  attack  mechanism  on  the  top  platform.  The  lift 
platform  is  supported  from  tlie  drag  platform  by  two  ver- 
tical plates  that  flex  only  in  the  lift  direction.  The  lift 
and  drag  platforms  are  also  connected  with  a flexure  with 
bonded  foil  strain  gauges  mounted  on  it.  The  drag  plat- 
form is  supported  by  two  vertical  plates  that  flex  only  in 
the  drag  direction  and  hang  from  two  more  vertical  flex- 
ible plates  attached  to  the  base  platform  of  the  balance. 
The  base  and  drag  platforms  are  also  connected  by  a flex- 
ure with  strain  gauges  mounted  on  it.  For  this  balance  a 
second  set  of  flexures,  for  both  lift  and  drag,  are  engaged 
when  the  loads  are  large.  For  the  range  of  forces  mea- 
sured in  this  investigation,  the  second  set  of  flexures  was 
never  engaged.  Figure  11  shows  a schematic  of  the  old 
balance  setup  in  the  wind  tunnel.  The  arrangement  with 
two  endplates  shown  in  Figure  11  is  known  as  arrange- 
ment number  1.  Arrangement  number  2,  not  shown,  has 
the  lower  endplate  removed  for  the  semi-infinite  tests. 

Thin-plate  models  for  current  investigation 

Keeping  in  mind  the  objective  of  this  first  phase 
of  tlie  investigation  which  was  to  study  the  aerodynamic 
characteristics  of  small,  low  aspect-ratio  flat  and  cam- 
bered wings,  several  thin,  flat  and  cambered  rectangu- 
lar aluminum  models  with  a thickness-to-chord  ratio  of 
1.93%  were  built.  Thin  models  were  selected  because 
birds  and  insects  have  very  tliin  wings.  The  models  either 
had  a 5-to-l  elliptical  leading  edge  and  a 3°  tapered  trail- 
ing edge,  or  5-to-l  elliptical  leading  and  trailing  edges. 
The  cambered  models  had  a circular  arc  shape  with  4% 
camber.  The  semi-span  aspect  ratios  (sAR)  tested  var- 
ied between  0.50  and  3.00.  The  root-chord  length  of  the 
models  was  either  4 in  (10.2  cm)  or  8 m (20.3  cm).  Fig- 
ure 12  shows  schematics  of  the  airfoil  geometries  for  the 
wings  with  a tapered  trailing  edge,  while  Table  1 gives 
the  dimensions  of  the  different  models  used.  With  the 
nomenclature  used  for  the  wing  designation,  the  first  four 
characters  define  the  nominal  dimensions  of  the  model. 
For  instance,  C8S4  means  a Chord  of  8 in  and  a Span  of 
4 in.  The  following  characters,  if  any,  define  the  shape: 
a C means  a cambered  plate  and  E means  an  elliptical 
trailing  edge  instead  of  a tapered  trailing  edge.  A max- 
imum span  of  12  m (30.5  cm)  was  chosen  so  that  these 
models  could  be  used  in  both  the  wind  and  water  tunnels. 


Figure  11:  UND-FB  1 balance  arrangement  (1)  with  two 
endplates  in  the  wind  tunnel 


thickness^  t 
(a)  Flat  plate 


(b)  Cambered  plate 


Figure  12:  Airfoil  geometry  for  models  with  tapered 
trailing  edge 
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Designation 

Chord  (in) 

Span  (in) 

sAR 

Thickness  (in) 

Camber  (%) 

C8S4 

7.973 

3.998 

0.5 

0.155 

0 

C8S8 

7.973 

8.003 

1.0 

0.154 

0 

C8S12 

7.985 

12.01 

1.5 

0.157 

0 

C4S8 

3.999 

8.019 

2.0 

0.077 

0 

C4S12 

4 

12.014 

3.0 

0.077 

0 

C8S4C 

7.975 

3.995 

0.5 

0.156 

4 

C8S8C 

7.983 

8 1 

1.0 

0.156 

4 

C8S12C 

7.908 

12.013 

1.5 

0.156 

4 

C4S8C 

3.995 

8 

2.0 

0.078 

4 

C4S12C 

3.936 

11.998 

3.0 

0.079 

4 

C8S12E 

7.969 

12.011 

1.5 

0.156 

0 

C8S12CE 

7.931 

12.011 

1.5 

0.157 

! 4 

C:  cambered;  E;  elliptical  trailing  edge 
Table  1 : Wi / ig  d imeus ions 


Tunnel  configurations 

Endplates  were  mounted  in  tlic  wind  and  water  tun- 
nels. The  plates  could  be  removed  to  simulate  citlier  a 
semi-infinite  model  or  a finite  model.  All  wings  tested 
were  held  at  the  quarter-chord  point  and  the  sting  was  cov- 
ered by  a streamlined  sting  covering  in  the  wind  tunnel 
and  a cylindrical  covering  in  the  water  tunnel.  The  gaps 
between  tlie  wing  and  tlie  endplates  were  adjusted  to  ap- 
proximately 0.03  m (0.8?7i/yv).  Mueller  and  Bums  (1982) 
showed  that  gap  sizes  varying  between  O.i  mm  and 
1.4  Trim  are  usually  acceptable  and  do  not  affect  the  re- 
sults. Furtliennore,  Rae  and  Pope  (1984)  suggest  that  the 
gap  be  less  than  0.005  x span.  For  a 12in  (30.5cm)  span 
model,  this  corresponds  to  a maximum  gap  size  of  0.06 zn 
(1.5  mm),  which  is  larger  than  tlie  gap  used  in  the  cur- 
rent investigation.  All  2D  tests  (or  infinite  wing/airfoil 
tests)  were  perfomied  with  both  endplates  present.  For 
semi -infinite  wings  (denoted  by  the  semi -span  aspect  ra- 
tio symbol  sAR),  the  bottom  plate  was  removed.  Finally, 
for  finite  wing  tests  (denoted  by  the  aspect  ratio  symbol 
A R),  both  endplates  and  the  sting  covering  were  removed. 

Data  acquisition  system  with  UND-FB1  balance 

Signals  from  tlie  strain  gauges  were  measured  with 
very  sensitive  instrumentation.  The  strain  gauges  were 
configured  in  a full  Wheatstone  bridge.  An  excitation 
voltage  of  5 V was  used  for  all  the  strain  gauge  bridges. 
The  bridge  signals  were  read  witli  an  instmmaitation  am- 
plifier circuit,  witli  available  gains  from  1 to  8,000.  The 
amplified  analog  signals  were  sent  to  the  computer  where 
they  were  then  converted  using  a four-chamiel,  12-bit  A/D 
converter  from  United  Electronic  Industries  (UEI).  Four 
data  chamicls  (lift,  drag,  moment  and  dynamic  pressure) 
could  be  measured.  All  tlie  data  was  acquired  using  a PC- 
based  data  acquisition  system  rumiing  tlie  LABVIEW® 
5 graphical  programming  language.  The  angle  of  attack 
was  controlled  manually  with  the  UND-FBl  balance. 


Procedure  for  data  acquisition 

Before  measuring  any  aerodynamic  force  and  mo- 
ment witli  eitlier  balance,  the  amplifier  gains  were  ad- 
justed to  maximize  the  output  signals  tliat  were  expected 
during  a given  set  of  experiments.  The  balance  was  then 
calibrated  using  known  masses.  The  lift,  drag  and  mo- 
ment axes  were  all  independent. 

For  tests  looking  at  the  aerodynamic  characteristics 
as  a function  of  angle  of  attack,  the  tunnel  velocity  was 
adjusted  with  the  model  at  a = 0°  to  yield  tlie  desired 
nominal  Reynolds  numba.  The  angle  of  attack  was  then, 
in  general,  set  to  a = —15°.  Data  was  taken  for  angles 
of  attack  up  to  a large  positive  angle  by  an  increment  of 
1°.  The  wing  was  dien  brought  back  to  a = 0°  by  an 
incremait  of  — 1°  in  order  to  see  if  hysteresis  was  present. 
Offset  readings  were  measured  for  all  four  data  acquisi- 
tion channels  before  the  tumiel  was  turned  on  with  the 
model  at  a = 0°.  At  the  end  of  the  nm,  the  tumiel  was 
turned  off  witli  tlie  model  at  a = 0°  and  drift  readings 
were  obtained  for  all  chamiels.  The  offset  voltage  for  a 
given  chaimel  was  subtracted  from  all  tlie  voltage  read- 
ings for  that  channel.  A percentage  of  the  drift  was  also 
subtracted  from  all  the  readings.  A linear  behavior  was 
assumed  for  the  drift.  This  means  that  if  n angles  of  at- 
tack were  tested  with  the  tunnel  running,  1/n  x di  ift  was 
subtracted  from  the  first  point,  2/n  x drift  was  subtracted 
from  tlie  second  point,  and  so  fortli.  Other  procedures  re- 
lated to  specific  applications  will  be  presented  in  the  text 
when  appropriate. 

Measurement  uncertainty  with  balance  UND-FB1 

Uncertainties  in  the  measurements  were  computed 
using  the  Kline-McClintock  tccluiiquc  (Kline  and  Mc- 
Clintock,  1953)  for  error  propagation.  The  two  main 
sources  of  uncertainty  were  the  quantization  error  and 
the  uncertainty  arising  from  the  standard  deviation  of 
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a given  mean  output  voltage.  The  quantization  error  is 
eg  = ^ volts  ^ jg  number  of  bits 

of  the  A/D  converter.  Optimizing  the  range  of  the  output 
voltages  can  help  to  reduce  the  uncertainties.  If  the  gain 
is  increased,  the  standard  deviation  of  the  mean  will  also 
be  increased,  but  tlie  ratio  of  the  standard  deviation  to 
the  mean  will  basically  remain  the  same.  However,  the 
uncertainty  from  the  quantization  error  will  be  reduced 
because  the  quantization  error  is  a fixed  value  (a  ftinction 
of  tlie  range  and  die  resolution  of  the  A/D  converter). 
The  ratio  of  the  quantization  error  to  the  mean  voltage 
will  dien  be  smaller  if  a larger  gain  is  used  and  a larger 
balance  output  mean  voltage  is  obtained. 

The  xmcertainty  in  die  angle  of  attack  was  deter- 
mined to  be  on  the  order  of  0.2°  ~ 0.3°.  Figures  13 
through  15  show  an  example  of  uncertainties  obtained  at 
Rec  = 60, 000  with  die  cambered  plates.  Error  bars  indi- 
cate the  uncertainty  in  Cl,  Cd  and  Cm/ 4^  The  average 
uncertainties  from  a ~ 3°  and  up  are  approximately  6% 
to  7%  for  Cl  and  Cp  and  10%  for  Cm/4- 


New  force/moment  aerodynamic  balance 
UND-FB2 


Description 

A new  platform  force/moment  balance  was  designed 
by  Matt  Fasano,  Professional  Specialist  at  the  Hessert 
Center  for  Aerospace  Research,  and  built  for  the  aerody- 
namic studies  on  low  aspect-ratio  wings  down  to  chord 
Reynolds  numbers  of  20,000.  The  design  of  this  new 
balance  (UND-FB2)  was  based  on  the  existing  balance 
(UND-FBl)  and  measures  lift,  drag,  and  pitching  moment 
about  die  vertical  axis.  It  is  an  external  balance  placed 
on  top  of  the  test  section  for  eidier  of  die  two  low-speed 
wind  tunnels  or  the  water  tunnel.  Due  to  the  better  sen- 
sitivity of  the  newly  designed  balance  (UND-FB2),  only 
tliis  balance  is  now  used  with  the  water  tunnel. 

With  this  balance,  lift  and  drag  forces  are  transmit- 
ted dirough  die  sting  which  is  mounted  directly  to  the 
moment  sensor  (see  Figure  16  for  a schematic  of  the  new 
balance).  The  moment  sensor  is  rigidly  mounted  to  the 
adjustable  angle  of  attack  mechanism  on  the  top  plat- 
form. The  lift  platform  is  supported  from  a platform, 
called  the  drag  platform,  by  two  vertical  plates  that  flex 
only  in  the  lift  direction.  The  lift  and  drag  platforms 
are  also  connected  with  a 1/8  in  x 3/8  in  x 1.5  in 
(3.2  7nm  X 9.5  mm  x 38.1  mm)  flexure  with  bonded 
foil  strain  gauges  mounted  on  it.  The  drag  platform  is 
supported  by  two  vertical  plates  that  flex  only  in  the 
drag  direction  and  hang  from  two  more  vertical  flexi- 
ble plates  attached  to  the  base  platform  of  the  balance. 
The  base  and  drag  platforms  are  also  connected  by  a 
l/16in  x 3/8in  x 1. Sin  (1.6mm  x 9.5mm  X 38.1mm) 
flexure  widi  strain  gauges  mounted  on  this  drag  flexure. 
Both  flexures  act  like  cantilever  beams  when  loads  are 
applied  to  the  balance. 


a (degrees) 

Figure  13:  Uncertainties  in  lift  coefficient  for  cambered 
plates  at  Rcc  = 60, 000  with  UND-FBl 


a (degrees) 

Figure  14:  Uncertainties  in  drag  coefficient  for  cambered 
plates  at  Rcc  ==  60,000  with  UND-FBl 


a (degrees) 

Figure  15:  Uncertainties  in  pitching  moment  coefficient 
for  cambered  plates  at  Rcc  = 60, 000  with  UND-FBl 
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Figure  16:  Schematic  of  the  new  balance  UND-FB2 


As  mentioned  earlier,  endplates  were  mounted  in 
botli  wind  and  water  tunnels.  Figure  17  shows  a 
schematic  of  the  new  balance  witli  the  endplates  in  place 
in  the  water  tunnel.  AH  wings  tested  were  held  at  the 
quarter-chord  point  and  the  sting  was  covered  by  a stream- 
lined sting  covering  in  the  wind  tunnel  and  a cylindrical 
covering  in  the  water  tunnel. 

The  moment  sensor  is  a Transducer  Techniques 
RTS-25  reaction  torque  sensor.  This  torque  sensor  uses 
bonded  foil  strain  gauges  and  is  rated  at  llyrnV /V  output. 
The  maximum  rated  capacity  is  25  oz  • cm  (J7.7  • cm) 

and  a torsional  stiffiicss  of  1, 32^1  N • cm/rad.  The  mo- 
mcnl  sensor  is  attached  to  an  adjustable  angle  of  attack 
mechanism  powered  by  a scr\^omotor  with  a controller. 

Electronics 

Signals  from  the  strain  gauges  are  measured  with 
very  sensitive  instrumentation.  The  strain  gauges  for  the 
drag  and  lift  flexures  arc  350  ohms  witli  a G-factor  of 
2.09  and  are  configured  in  a fiill  Wheatstone  bridge.  An 
excitation  voltage  of  5 V is  used  for  all  the  strain  gauge 
bridges.  The  bridge  signals  are  read  with  an  instrumenta- 
tion amplifier  circuit  witli  a gain  as  higli  as  8,000.  Due  to 
the  sensitivity  of  the  circuit  many  precautions  were  made 
to  reduce  noise.  At  first,  Ni-Cd  rechargeable  batteries 
were  used  to  power  the  amphfiers,  analog-to-digital  con- 
verters, and  the  excitation  voltage  for  the  strain  gauges. 
A DC  power  supply  is  now  being  used  because  of  the 
quick  discharge  of  the  batteries  during  data  acquisition. 


Figure  17:  New-halance  arrangement  (I)  in  the  water 
tunnel  with  two  endplates 


The  batteries  used  could  not  provide  a constant  voltage 
for  several  lioxus. 

The  amplifiers  and  analog-to-digital  converters  are 
mounted  on  a circuit  board  placed  in  a control  box  with 
switches  and  potentiometers  (pots)  to  adjust  tlie  gains, 
offsets  and  balance  tlie  Wheatstone  bridges.  Four  data 
channels  (lift,  drag,  moment  and  dynamic  pressure  when 
necessary)  can  be  measured  quasi-simultaneously.  The  in- 
put differaitial  signal  from  each  channel  is  sent  through 
two  amplifiers  from  Analog  Devices:  a precision  instm- 
meiitation  AD624  amplifier  (gain  of  I,  100,  200  or  500) 
and  a software  programmable  AD526  gain  amplifier  (gain 
of  1,  2,  4,  8 or  16).  The  amplified  analog  single-ended 
signals  from  tlie  four  channels  are  then  converted  to  dig- 
ital signals  using  a Burr- Brown  AD7825,  four-chamiel, 
16-bit  analog-to-digital  converter.  The  signals  are  then 
salt  to  a National  Instruments  digital  data  acquisition  card 
(NIDAQ)  in  a data  acquisition  computer.  The  amplified 
single-ended  analog  signals  can  also  be  sent  directly  to 
the  computer,  thus  bypassing  tlie  16-bit  A/D  converters; 
the  signals  are  tlien  converted  using  tlie  four-chaimel,  12- 
bit  A/D  UEI  converter,  mentioned  earlier.  With  tlie  16-bit 
A/D  system,  each  amplifier  circuit  is  identical  except  for 
the  fourtli  chamiel  where  the  amplifiers  can  be  bypassed 
and  a single-ended  signal  can  be  sent  directly  to  the  16- 
bit  A/D  converter.  Figure  18  is  a simplified  schematic  of 
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the  electronic  circuitry  used. 


Offset 


Figure  18:  UND-FB2  balance  electronics 


Data  acquisition 

All  tlie  data  was  acquired  using  a PC-based  data 
acquisition  system  running  the  LABVIEW®  5 graphical 
programming  language.  The  NIDAQ  card  was  first  used 
for  data  acquisition.  The  UEI  card  is  now  used  with  the 
new  balance  when  severe  noise  interferes  with  the  data 
and  the  NIDAQ  card  cannot  be  used. 

The  data  acquisition  process  with  the  new  balance 
was  automated.  The  angle  of  attack  can  be  automatically 
varied  from  a pre-determined  list  of  angles  of  attack.  The 
range  is  usually  adjusted  in  order  to  be  able  to  observe 
stall. 

Specifications 

Since  the  force/moment  balance  includes  very  sensi- 
tive flexures  and  strain  gauges,  the  applied  forces  and  mo- 
ment cannot  exceed  certain  limits.  These  limiting  forces 
and  moment  were  determined  conservatively  and  are  listed 
in  Table  2.  The  limiting  forces  and  moment  for  the  UND- 
FBl  balance  are  also  included  for  comparison  (see  Huber, 
1985).  In  order  to  be  able  to  move  the  balance  and  mount 
the  models  without  permanently  deforming  the  flexures, 
locking  pins  are  used  to  restrain  the  balance.  These  lock- 
ing pins  must  be  removed  when  taking  data. 

Sources  of  noise 

For  all  measurements,  digital  filtering  in 
LABVIEW®  was  necessary  to  reduce  noise  gener- 
ated by  the  servomotor  used  to  change  the  angle  of 


TJND-FB2 

UND-FBl 

Positive  lift 

7.0 

39.2  7V 

Negative  lift 

-3.5 

-14.7  7V 

Positive  drag 

2.0  7V 

14.7  7V 

Moment 

1 15.0  N • cm 

1 226  N • cm 

Table  2:  Maximum  force/moment  balance  specifications 


attack,  and  also  the  motor  of  tlie  water  tunnel  when  water 
tunnel  tests  were  performed.  A low-pass  Butterworth 
filter  with  a cut-off  frequency  of  50  Hz  and  of  order 
5 was  used.  A study  on  the  effect  of  the  filter  and  its 
order  showed  tliat  the  mean  voltages  were  basically  not 
affected,  but  the  standard  deviations  of  the  means  were 
greatly  reduced.  It  was  discovered  during  preliminary 
calibrations  with  the  NIDAQ  card  that  the  servomotor 
was  causing  noise  in  die  data.  With  the  motor  ON,  the 
standard  deviations  of  the  samples  (4,000  data  points 
measured  at  a sampling  frequency  of  500  Hz)  were 
larger  than  those  without  the  motor  ON,  although  the 
mean  values,  dius  the  calibration  coefficients,  were  the 
same.  It  was  found  that  isolating  the  motor  from  the 
balance  helped  to  reduce  the  standard  deviations.  A diin 
plastic  sheet  was  then  placed  between  the  aluminum 
motor  support  and  the  aluminum  top  plate,  lift  platform, 
of  the  balance.  Moreover,  plastic  screws  were  used  to 
mount  the  motor  support  to  the  balance.  This  eliminated 
any  aluminum/aluminum  contact  between  the  motor  and 
the  balance.  Noise  caused  by  a metal-to-metal  contact 
between  a force/moment  balance  and  a motor  was  also 
detected  in  a previous  investigation  at  the  University 
of  Notre  Dame  (Pelletier,  1998).  Figure  19  shows  the 
standard  deviations  of  the  samples  for  a lift  chamiel 
calibration  example  with  and  without  the  isolation 
plastic. 


'o 


0.5 


0.4 


0.3 


o With  motor  OFF 

^ With  motor  ON  (no  isolation) 
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Figure  19:  Effect  of  the  motor  on  the  standard  deviation 
with  NIDAQ  data  acquisition  card 
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Calibration 


After  construction  of  the  balance  and  electronics, 
calibrations  were  performed  to  ascertain  the  linearity  of 
the  balance  for  all  its  independent  axes:  lift,  drag  and 
moment.  These  calibrations  were  perfomicd  by  moving 
the  upper  plate,  lift  platform,  of  the  balance  or  by  apply- 
ing  a torque  to  tlie  moment  sensor  by  placing  precision 
weights  of  known  mass  in  a container  connected  to  the 
plate,  or  moment  sensor.  The  container  was  comiected  to 
the  upper  plate  by  running  a string  over  the  calibration 
pulley  aligned  with  the  axis  to  calibrate  (for  lift  and  drag) 
or  perpendicular  to  tlic  torque  applicator  (pin  comiected 
to  the  sting  used  to  apply  a known  torque)  for  tlie  mo- 
ment calibration.  For  tliis  moment  calibration,  the  lift  and 
drag  locking  pins  were  used  to  prevent  movement  in  the 
lift  and  drag  direction  as  a moment  was  applied  to  the 
sensor.  Several  calibrations  were  perfonned  to  look  for 
repeatability  and  linearity.  Figures  20  tlirougli  22  show 
examples  of  tlie  calibration  curves  tliat  were  repeatedly 
obtained. 
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Figure  20:  Lift  calibration  for  the  new  balance 
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UND-FB2  performance 

Once  die  new  balance  and  its  electronics  were  built 
and  calibrations  had  bcai  perfonned,  it  was  tested  to  see 
if  tlie  results  compared  to  published  data.  All  results  pre- 
sented in  Uiis  paper  have  beai  corrected  for  solid  block- 
age, wake  blockage  and  streamlined  curv^ature  using  tech- 
niques presented  by  Pankliurst  and  Holder  (1952)  and  Rae 
and  Pope  (1984),  A series  of  two-dimensional  tests  were 
conducted  on  different  models. 

Models 

In  die  balance  validation  phase,  two  circular  cylin- 
ders were  tested.  The  diameters  of  the  two  cylinders  were 
0.75  m (1.9  cm)  and  1.255  in  (3.2  cm)  and  they  both  had 
a length  of  12  in  (30.5  cm).  An  Eppler  61  airfoil  model, 
whose  profile  was  shown  in  Figure  4,  was  also  used  to 
test  the  balance.  The  model  also  had  a length  of  12  in 
(30.5  cm)  and  a chord  of  4.906  in  (12.5  cm). 

Cylinder  results 

The  new  balance  was  first  tested  by  measuring  the 
two-dimensional  drag  on  two  circular  cylinders  in  the 
low-speed  wind  tunnel.  Figure  23  shows  die  drag  co- 
efficient as  a function  of  Reynolds  number.  Results  of 
the  present  investigation  were  compared  to  results  by 
Wieselsberger,  digitized  from  Boundaiy-Layer  Theory  by 
Schlichting  (1979).  There  is  a good  agreement  between 
the  two  sets  of  data. 

Eppler  61  airfoil  results 

The  balance  was  then  tested  by  measuring  die  two- 
dimensional  lift  and  drag  on  die  Eppler  61  airfoil.  Results 


Figure  21:  Drag  calibration  for  the  new  balance 


Calibration  voltage  (volts) 

Figure  22:  Pitching  moment  calibration  for  the 
new  balance 
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Re 

Figure  23:  2D  drag  coefficient  of  two  circular  cylinders 
with  UND-FB2 


were  obtained  in  the  wind  tunnel  and  the  water  tunnel  for 
several  Reynolds  numbers.  Figures  24  and  25  show  the 
two-dimensional  lift  and  drag  coefficients  respectively  in 
the  wind  tunnel  for  different  nominal  Reynolds  numbers, 
i.e.,  values  used  to  adjust  the  velocity  in  the  tunnel  with 
the  model  at  a = 0°. 

Results  for  Q indicate  a significant  difference  be- 
tween results  at  Rcc  > 60,000  and  those  for  Rcc  < 
60, 000.  For  large  Re^.,  Ci  increases  smoothly  with  angle 
of  attack  a.  For  smaller  Re,,,  the  lift-curve  slope  Ci^ 
is  smaller  for  2°  < a < 8""  and  there  is  a sharp  rise 
in  Cl  at  a 8°.  Similar  results  have  been  obtained  by 
Altliaus  (1980)  and  shown  in  Figure  26.  Althaus  used 
a strain  gauge  balance  arrangement  to  measure  lift  and  a 
wake  rake  to  measure  drag.  A drawback  of  using  a wake 
rake  will  be  addressed  later.  Althaus  did  observe  a small 
hysteresis  loop  at  low  Reynolds  numbers.  No  apparent 
hysteresis  was  observed  in  tlie  current  study. 

The  sharp  rise  in  Ci  at  low  Reynolds  numbers  is 
believed  to  be  tlie  result  of  a laminar  separation  bub- 
ble on  tlie  upper  surface  of  the  wing.  O’Meara  and 
Mueller  (1987)  showed  that  the  length  of  the  separation 
bubble  tends  to  increase  with  a reduction  in  Re^.  A re- 
duction in  the  turbulence  intensity  also  tends  to  increase 
the  length  of  the  bubble.  The  lift-curve  slope  is  affected 
by  separation  bubbles.  A longer  bubble  is  usually  associ- 
ated widi  a decrease  in  tlie  lift-curve  slope  (Bastedo  and 
Mueller,  1985).  This  is  the  kind  of  behavior  observed 
witli  tlie  Eppler  61  airfoil  in  this  investigation.  Flow 
visualization  by  Mueller  and  Bums  (1982)  showed  the 
presence  of  a separation  bubble  on  the  Eppler  61  airfoil 
at  Rcc  46, 000. 

In  the  water  tunnel,  the  turbulence  intensity  is  larger 
than  in  the  wind  tunnel.  Therefore,  the  smaller  Ci^  and 
sharp  rise  in  Ci  observed  in  tlie  wind  tunnel  for  2°  < 
Q < 8°  might  not  be  present  at  all  in  the  water  tunnel 
data  for  the  same  Reynolds  numbers.  This  is  exactly 


a (degrees) 

Figure  24:  2D  lift  coefficient  on  the  Eppler  61  airfoil  in 
the  wind  tunnel  with  UND-FB2 


CL  (degrees) 

Figure  25:  2D  drag  coefficient  on  the  Eppler  61  airfoil 
in  the  wind  tunnel  with  UND'FB2 


what  happened  for  Re^  = 42, 000  and  Rcc  — 46, 000,  as 
shown  in  Figure  27.  The  drag  coefficient  appears  to  be 
less  affected,  within  the  uncertainty  of  the  measurements. 
Figure  28  shows  the  Cd  curves  obtained  in  the  water 
tunnel  for  the  Eppler  61  airfoil. 

Figures  29  through  31  show  comparisons  of  the  cur- 
rent Eppler  61  results  with  published  data.  There  is,  in 
general,  a good  agreement  betweai  the  current  data  and 
published  data.  The  most  significant  difference  is  in  the 
stall  angle;  there  appears  to  be  a 2°  difference  in  astaii- 

Measurement  uncertainty  with  UND-FB2 

Uncertainties  in  the  measurements  were  computed 
using  the  Kline-McClintock  technique  (Kline  and  Mc- 
Clintock,  1953)  for  error  propagation.  As  indicated  ear- 
lier, tlie  two  main  sources  of  uncertainty  were  tlie  quanti- 
zation error  and  the  uncertainty  arising  from  tlie  standard 
deviation  of  a given  mean  output  voltage.  The  quantiza- 
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Figure  26:  Althaus'  results  for  2D  lift  coefficient  on 
the  Eppler  61  airfoil  (Althaus,  1980) 


tion  error,  described  earlier,  was  smaller  witli  the  NIDAQ 
card  lhaii  willi  the  UEI  card  due  to  die  better  resolution 
of  the  A/D  converter  (16  bits  compared  to  12  bits).  The 
uncertainty  in  the  angle  of  attack  was  determined  to  be  on 
the  order  of  0.2°  ~ 0.3°.  The  error  from  the  encoder  was 
negligible.  The  encoder  offered  an  excellent  resolution 
of  2,000  counts  per  degree,  which  gave  an  uncertainty  of 
5 X lO”"*.  Figures  32  and  33  show  a comparison  of  wind 
tunnel  and  water  tunnel  results  for  tlie  Eppler  61  airfoil 
at  Rcc  = 42, 000.  Error  bars  indicate  the  uncertainty  in 
Cl  and  Cd  when  the  NIDAQ  card  is  used.  The  average 
uncertainties  for  Ci  and  Cd  in  tlie  range  of  angles  of  at- 
tack tested  arc  approximately  4%  in  tlie  water  tumiel  and 
6%  in  die  wind  tumiel. 

The  new  balance  in  itself  is  also  more  sensitive  dian 
the  old  balance.  This  allows  experiments  at  smaller  ve- 
locities in  the  water  tumiel.  As  of  now,  results  with  dif- 
ferent wings  have  shown  a high  degree  of  repeatability  for 
Reynolds  numbers  as  low  as  40, 000.  A major  challenge 
in  measurements  at  Reynolds  numbers  below  40,000  is 
being  able  to  measure  drag  accurately.  At  Rcr  = 20, 000, 
the  minimum  drag  can  be  as  low  as  0.02  A,  which  corre- 
sponds to  a load  of  approximately  2 grams.  A fine  drag 
calibration  of  die  balance  showed  that  1 gram  was  suf- 
ficient to  deflect  the  drag  flexure  and  yield  a reasonable 
output  voltage.  However,  diis  deflection  is  often  on  the 
order  of  sigial  noise. 
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Figure  27:  2D  lift  coefficient  on  the  Eppler  61  airfoil  in 
the  water  tunnel  with  UND-FB2 
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Figure  28:  2D  drag  coefficient  on  the  Eppler  61  airfoil 
in  the  water  tunnel  with  UND-FB2 
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Results  for  thin  wings 

This  section  will  present  results  for  thin  flat  and 
cambered  wings.  Some  additional  issues  associated  with 
detennining  aerodynamic  characteristics  as  a function  of 
Reynolds  numbers  will  also  be  addressed.  Accurate  mea- 
surements of  Cl  and  Cd  with  endplatcs  and  small  aspect- 
ratio  models  are  difficult  to  obtain  at  low  Reynolds  num- 
bers because  of  die  interaction  between  die  thick  bound- 
ary layers  on  die  cndplates  and  die  flow  around  the  wing, 
which  results  into  a tliree-dimensional  flow  along  the  span 
of  die  model.  This  must  be  kept  in  mind  when  examining 
the  following  results. 

Flat-plate  wings 

Some  results  for  the  flat-plate  models  (two  endplates 
for  2D  tests  and  one  endplate  for  semi  aspect-ratio  tests) 
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Figure  29:  Comparison  of  2D  lift  coefficient  on  the 
Eppler  61  airfoil  and  published  data  with  UND-FB2 


Re=  42,000 


Figure  32:  Comparison  of  the  2D  lift  coefficient  on  the 
Eppler  61  airfoil  in  wind  and  water  tunnels  with  UND-FB2 
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Figure  30:  Comparison  of  2D  drag  coefficient  on  the 
Eppler  61  airfoil  and  published  data  with  UND'FB2 
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Figure  33:  Comparison  of  the  2D  drag  coefficient  on  the 
Eppler  61  airfoil  in  wind  and  water  tunnels  with  UND'FB2 
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Figure  31:  Comparison  of  2D  drag  polar  for  the  Eppler 
61  airfoil  and  published  data  with  UND-FB2 


can  be  seen  in  Figures  34  through  39  for  Re^  = 80, 000 
and  Rcc  = 140, 000.  Figures  34  and  37  show  a signifi- 
cant reduction  in  the  lift-curve  slope  for  semi-infinite 
wings. 

The  lift-curve  slope  values  obtained  from  the  wind 
tunnel  data  are  compared  to  theoretical  values  for  thin 
wings  of  different  semi-span  aspect  ratios  in  Figure  40. 
Equation  1 from  Anderson  (1991)  was  used  to  estimate 
the  theoretical  values  of  : 


oo 


(1) 


where  ao  is  the  2D  lift-curve  slope  in  1/degrees,  AR  is 
the  aspect  ratio  of  the  full  wing  (AR  = 2 * sAR)  and  r 
is  the  Glauert  parameter  (equivalent  to  an  induced  drag 
factor)  varying  typically  between  0.05  and  0.25.  The  2D 
value  ao  was  determined  to  be  ao  = 0.0938/deg.  This 
corresponds  to  the  average  of  all  the  slopes  (for  all 
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Figure  34:  Lift  coefficient  on  flat  plates  at 
- 80,  000  with  UND-FBl 
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Figure  35:  Drag  coefficient  on  flat  plates  at 
Re,,  - 80,  000  with  UND^FB  I 
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Figure  37:  Lift  coefficient  on  flat  plates  at 
Re„  = 140,000  with  UND-FBl 
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Figure  38:  Drag  coefficient  on  flat  plates  at 
Re„  = 140, 000  with  UND-FBl 
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Figure  36:  Pitching  moment  coefficient  on  flat  plates 
at  Re,,  - 80, 000  with  ' UND-FBl 
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Figure  39:  Pitching  moment  coefficient  on  flat  plates  at 
Rcc  - 140, 000  with  UND-FBl 
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Figure  40:  Lift-curve  slope  for  flat-plate  models  in  wind 
tunnel  with  UND-FBl 


Reynolds  numbers  considered)  for  an  infinite  aspect  ra- 
tio (1/sAR  ^ 0).  This  value  was  picked  instead  of  the 
conventional  value  of  ao  — 27r/rad  = 0.1/deg  given  by 
thin-airfoil  theory.  Figure  40  shows  a very  good  agree- 
ment between  the  experimental  values  of  and  the 
theoretical  values  estimated  by  Equation  1. 

As  tlie  aspect  ratio  was  decreased,  Figures  34  and  37 
also  show  tliat  the  linear  region  of  tlie  vs  a curve 
became  longer  and  agtaii  tended  to  increase.  Moreover, 
both  figures  show  that  there  was  no  abrupt  stall  for  low 
aspect-ratio  wings.  For  these  low  aspect  ratios,  Cl  often 
reached  a plateau  and  then  remained  relatively  constant, 
or  even  started  to  increase,  for  increasing  angles  of  attack. 

Changing  the  aspect  ratio  of  the  models  did  not 
appear  to  have  a measurable  effect  on  the  drag  coef- 
ficient at  Rcc  = 80,000,  as  shown  in  Figure  35.  At 
Rcc  = 140, 000,  increasing  the  aspect  ratio  had  the  unex- 
pected effect  of  increasing  Co  for  angles  greater  than  5°. 
No  measmable  differoice  was  encountered  in  the  range 
-5°  < a < 5°. 

Finally,  Figures  36  and  39  show  the  pitching  mo- 
ment at  the  quarter  chord.  Both  figures  indicate  a slightly 
positive  slope  Cm^  around  a = 0°,  even  when  consider- 
ing the  uncertainty.  This  would  imply  that  the  flat-plate 
models  were  statically  unstable  around  a = O'^ . Increas- 
ing the  Reynolds  number  from  80, 000  to  140, 000  tended 
to  reduce  the  slope  of  C^/4.  The  model  with  a semi- 
span aspect  ratio  of  3 indicated  an  irregular  behavior  at 
Rcc  = 140, 000  for  C^/4;  the  pitching  moment  was  not 
zero  at  a — 0°.  This  case  will  have  to  be  repeated 

Aerodynamic  characteristics  as  a function  of 
Reynolds  number:  a different  method 

For  tests  without  endplates,  another  balance  arrange- 
ment, denoted  arrangement  number  3,  was  used  and  is 
presented  in  Figure  41.  The  lift  and  drag  forces  measured 
by  the  balance  were  for  the  wing-sting  combination.  The 
lift  on  the  sting  was  basically  zero.  However,  the  drag 


of  the  sting  alone,  which  dominated  the  total  wing-sting 
drag,  was  not  zero  and  was  subtracted  from  the  wing-sting 
values  to  get  the  Cd  of  tlie  flat-plate  wing  alone. 


Figure  41:  Balance  arrangement  (3)  for  finite  wing  tests 
with  the  new  balance 


In  general,  when  investigators  try  to  determine  how 
Cl, and  vary  with  Reynolds  numbers,  they  deter- 
mine Cl„^^^  and  CL)„,i„from  Cl  vs  a and  Cp  vs  a curves 
at  different  Reynolds  numbers.  It  has  been  observed  in 
this  investigation  that  the  values  obtained  do  not  always 
match  the  expected  trend  for  drag  because  of  the  difficulty 
involved  in  measuring  the  very  small  drag  forces.  A slight 
offset  in  one  Cl  vs  a or  Cl>  vs  a curve  can  lead  to  jagged 
or  vs  Rec  curves.  A better  tech- 

nique was  found  to  obtain  Cl>„^j„  vs  Rec  (the  values  of 
CLm^x  Rec  are  of  lesser  importance  because  micro-air 
vehicles  will  rarely  fly  at  For  this  technique,  the 

angle  of  attack  was  fixed  to  the  angle  yielding  the  lowest 
Cd  in  a Cl)  vs  a curve,  and  measurements  were  taken  for 
a series  of  increasing  and  decreasing  Reynolds  numba*s 
without  stopping  the  tunnel.  Results  obtained  with  the 
new  balance  UND-FB2  using  this  technique  on  a finite 
wing  of  aspect  ratio  AR  ==  1 in  the  wind  tunnel,  pre- 
sented in  Figures  42  and  43,  are  promising  and  the  trends 
obtained  matched  the  expected  reduction  in  CD^^.^with 
increasing  Reynolds  numbers. 
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Figure  42:  Lift  coefficient  variation  with  Re  for 
AR  = 1 flat-plate  wing  with  UND-FB2  (C8S8) 
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(b)  a = 10° 

Figure  43:  Drag  coefficient  variation  with  Re  for 
AR  = I flat-plate  wing  with  UND-FB2  (C8S8) 


Results  at  tlie  angle  of  attack  for  ~ 3")  also 

show  an  increase  in  Ci  and  a decrease  in  Cp  witli  in- 
creasing Reynolds  numbo’.  The  results  of  Figure  42  at 
a = 10°  are  smaller  than  the  lift  coefficient  presented  in 
Figure  34  for  a model  witli  sAR  ~ 1 because  of  tlie  lack 
of  an  endplate.  For  a given  model,  adding  an  endplate 
leads  to  an  increase  in  lift  compared  to  the  case  without 
an  endplate,  as  shown  in  Figure  44.  Adding  one  end- 
plate did  not  have  a significant  effect  on  Cp,  as  shown 
in  Figure  45.  The  wing  used  for  this  series  of  tests  had  a 
nominal  chord  c ~ Sin  and  span  6 = 12 in  [CSS  12]. 
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Figure  44:  Lift  coefficient  variation  with  Re  (C8SI2) 
\>ith  UND-FB2  at  a = 3° 
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Figure  45:  Drag  coefficient  variation  with  Re  (CBS 1 2) 
with  UND-FB2  a/  a - 3° 

At  an  angle  of  attack  a ^ 0°,  adding  one  or  two 
endplates  did  not  affect  Ci  for  Re^  > 60,000;  it  re- 
mained basically  around  the  expected  value  of  zero,  as 
shown  in  Figure  46.  Moreover,  adding  endplates  did  not, 
once  again,  have  a measurable  effect  on  Cp,  as  presented 
in  Figure  47.  Since  Figure  47  is  for  a 0°,  it  represents 
versus  Reynolds  number  for  flat-plate  wings.  Ex- 
perimental results  are  also  compared  to  tlieory  (Blasius: 
" 2.66/?.e“2)  and  CFD  results  in  Figure  47.  The 
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CFD  results  were  computed  by  Greg  Brooks  {Air  Force 
Research  Laboratory,  Wright-Patterson  Air  Force  Base) 
using  COBALT,  a parallel,  implicit,  unstructured,  finite 
volume  laminar  CFD  code  based  on  Godunov’s  exact  Rie- 
mann  method,  developed  by  the  Air  Force  Research  Lab- 
oratory (see  Strang,  Tomaro  and  Grismer,  1999).  All  sets 
of  data  indicate  the  same  trend  The  e^erimental  data 
was  always  larger  than  theoiy  and  the  CFD  results.  This 
could  have  been  caused  by  surface  rou^mess,  imperfect 
flow  conditions,  and  so  forth.  Since  all  wind  tunnel  tests 
with  the  flat  and  cambered  wings  are  usually  conducted 
at  Reynolds  numbers  greater  than  60, 000,  the  results  pre- 
sented in  Figures  42  through  47  for  Re  < 60, 000  should 
be  analyzed  with  caution.  The  velocity  in  the  wind  tun- 
nel is  usually  too  low  for  i^e  < 60, 000  to  yield  reliable 
results. 
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Figure  46:  Lift  coefficient  variation  with  Re  (C8S12) 
with  UND-FB2  a/  a -=  0° 
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Figure  47:  Drag  coefficient  variation  with  Re  (CSS  12) 
with  IJND-FB2  a = 0° 


The  results  of  Figure  47  seem  to  indicate  that  for 
thin  flat-plate  wings  at  o;  = 0°  the  drag  coefficient  acting 
on  the  wing  is  independent  of  aspect  ratio.  Results  of 
Figure  47  were  then  compared  to  the  drag  coefficients 


Figure  48:  Minimum  drag  coefficient  with  Re  for 
AR  > 1 flat-plate  wings  with  UND-FB2 


for  the  AR  = 1 plate  and  no  difference  was  obtained, 
as  shown  in  Figure  48.  For  aspect  ratios  greater  than 
one,  the  minimum  drag  coefficient  is  then  independent 
of  model  size  and  the  presence  of  endplates.  For  lower 
aspect  ratios,  preliminary  tests  indicated  a larger 
More  work  is  in  progress  to  fully  understand  the  behavior 
of  C^^^.^^as  a function  of  Reynolds  numbers  for  the  very 
low  aspect-ratio  wings  (AR  < 1).  Similar  tests  will  also 
have  to  be  conducted  on  the  cambered  wings.  The  non- 
measurable  difference  in  Reynolds  numbers 

might  just  be  valid  for  flat-plate  wings.  Adding  camber 
might  change  the  results.  Preliminary  results  seem  to 
indicate  this  trend. 

Cambered-plate  wings 

Results  were  also  obtained  for  cambered-plate  mod- 
els using  the  balance  arrangement  with  one  or  two  end- 
plates (semi-infinite  or  2D  tests).  In  general,  camber 
led  to  better  aerodynamic  characteristics  due  to  an  in- 
crease in  lift,  even  though  drag  also  increased.  Figures  49 
through  54  show  some  results  for  the  cambered  plates  at 
Rce  = 60,000  and  Re^  = 140,000.  With  cambered 
plates,  C'o^.i^.was  slightly  larger  than  for  flat  plates.  The 
maximum  lift  coefficient  was  also  larger,  as  expected. 
Moreover,  the  variation  in  Cl  with  angle  of  attack  at 
small  angles  was  less  linear  for  cambered  plates  than  for 
flat  plates.  Finally,  the  behavior  of  the  moment  coefficient 
Cm/ 4 for  the  cambered  plates  was  very  differmt  than  the 
behavior  with  the  flat  plates.  A rise  in  Cm./ 4 occurred 
after  a ^ 5°,  leading  to  a hump  at  around  a ^ 10°.  This 
was  not  observed  with  the  flat  plates.  Flow  visualization 
will  hopefiiUy  explain  this  behavior. 

Equation  1 was  also  used  to  compare  the  experi- 
mental values  of  Cl^  at  aci^=o  for  the  cambered  plates 
to  theoretical  values.  The  2D  value  ao  used  was  ao  == 
0.1097/deg.  Figure  55  shows  a good  agreement  between 
theory  and  e:?q)eriments. 
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Figure  49:  Lift  coefficient  on  cambered  plates  at 
Rcr  = 60,  000  with  UND-FBl 
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Figure  52:  Lift  coefficient  on  cambered  plates  at 
Rcc  = 140',000  with  UND-FBI 
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Figure  50:  Drag  coefficient  on  cambered  plates  at 
Rcr.  - 60,  000  with  UND-FBl 
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Figure  53:  Drag  coefficient  on  cambered  plates  at 
Rcr^  = 140, 000  with  UND-FBl 
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Figure  51:  Pitching  moment  coefficient  on  cambered 
plates  at  ~ 60,000  with  UND-FBl 
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Figure  54:  Pitching  moment  coefficient  on  cambered 
plates  at  Re^^  = 140, 000  with  UND-FBl 
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Figure  55:  Lift-curve  slope  for  cambered-plate  models 
in  wind  tunnel  with  UND-FBl 


Since  the  cambered  wings  showed  better  aerody- 
namic characteristics,  and  hence  are  more  suitable  in  the 
design  of  micro-air  vehicles,  only  performance  data  for 
the  cambaed  plates  is  presented.  Figures  56  through  59 
show  the  behavior  of  Cr  , Cn  . , (^)  and 

I 1 as  a function  of  Reynolds  numb  a.  The  max- 

\ / max 

imum  L/D  ratio  is  related  to  the  maximum  range  for  a 
propeller  driven  airplane,  while  the  maximum  f Co 

is  related  to  best  endurance  (longest  flying  time  possi- 
ble). As  expected,  increased  with  Reynolds  number 

and  aspect  ratio  in  tlie  range  of  Reynolds  numbers  tested. 
The  same  expected  behavior  was  obtained  for  (;§)n^a^^and 

( other  hand,  showed  an  in- 

\ / max 

crease  with  decreasing  Reynolds  number,  as  was  also  ex- 
pected. The  maximum  L/D  generally  occurred  at  a ==  3° 

to  4°,  while  f ) occurred  at  a = 4°  to  5°.  It 

\ / max 

is  important  to  remember  that  the  endplates  have  been 
shown  to  have  an  effect  on  the  lift  coefficients  of  flat- 
plate  wings,  and  probably  cambered-plate  wings  also.  As 
mentioned  earlier,  the  effect  of  the  endplates  on  the  drag 
characteristics  of  cambered-plate  wings  is  still  under  in- 
vestigation. The  results  presented  in  this  report  must  then 
be  analyzed  with  the  possible  effect  of  tlae  endplates  in 
mind;  the  numerical  values  should  not  be  taken  as  the  ul- 
timate results.  The  effect  of  the  endplates  on  the  pitching 
moment  has  not  been  investigated  for  both  flat-plate  and 
cambered-plate  wings. 


Effect  of  trailing-edge  geometry 


Four  models  [c  = Sin  (c  ^ 20.3cm)  and  6 — 12m 
(b  — 30.5  cm)]  were  tested  in  the  wind  tunnel  at  sev- 
eral chord  Reynolds  numbers  to  see  if  the  trailing-edge 
geometry  had  any  influence  on  the  aerodynamic  charac- 
teristics of  flat  plates  and  cambered  plates  at  low  chord 


Re, 

Figure  56:  Maximum  lift  coefficient  as  a function  of 
Re c for  cambered  wings  with  UND-FBl 
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Figure  57:  Minimum  drag  coefficient  as  a function 
of  Rccfor  cambered  wings  with  UND-FBl 


Figure  58:  Maximum  L/D  ratio  as  a function  of  Re^ 
for  cambered  wings  with  UND-FBl 


8-23 


Figure  59:  Mcaimum  ICq  ratio  as  a fimction 
of  Rcr  for  cambered  wings  with  UND-FBl 


a (degrees) 

Figure  60:  Trailing'edge  geometry  effect  on  lift 
coefficient  at  = 80,000  on  flat  plates  in  the 
wind  tunnel  with  UND-FBl 


Reynolds  numbers.  The  first  two  models  had  a tapered 
trailing  edge,  while  tlie  otlier  two  models  had  an  ellipti- 
cal trailing  edge.  Results  were  obtained  for  infinite  mod- 
els (2D  case)  and  models  witli  a semi-span  aspect  ratio 
sAR  1.5.  For  botli  cases,  no  significant  difference  was 
observed  in  Ci  or  C/,  and  Cd  or  Ca,  as  a function  of 
trailing-edge  geometry,  as  shown  in  Figures  60  and  61  for 
Re,,  “ 80,000.  A differaice  was  however  observed  in 
the  moment  coefficient  For  a sharp  trailing  edge, 

Cfn„  often  appeared  to  be  positive  around  q = O'",  even 
with  the  uncertainty  considered  (error  bars  in  C^n/4 
about  the  size  of  the  symbols).  Witli  the  elliptical  trailing 
edge,  tlie  2D  cases  at  Re,.  = 80,000  showed  a stable 
negative  value  of  shown  in  Figure  62.  For  a 

semi-span  aspect  ratio  of  1.5,  was  basically  zero  at 
a = 0”.  Flow  visualization  to  be  perfonned  later  may 
explain  this  phenomenon. 

With  the  cambered  plates,  tliere  was  basically  no 
difference  between  a sharp  trailing  edge  and  an  elliptical 
trailing  edge  at  Re„  = 80,000,  as  shown  in  Figures  63 
through  65.  Results  with  the  cambered  plates  seem  to 
agree  with  Laitone  (Laitone,  1996  and  1997),  who  showed 
that  at  low  Reynolds  numbers,  a sharp  trailing  edge  is  not 
as  critical  as  for  larger  Reynolds  numbers. 

The  influence  of  the  leading-edge  geometry  was  also 
investigated  by  looking  at  tlic  lift  and  drag  characteristics 
of  a flat-plate  model  in  2D  and  sAR  = 1.5  configurations 
in  the  water  tunnel  at  Re,,  —■  39, 000  and  Re„  = 60, 000. 
For  these  tests,  the  existing  CSS  1 2 model  was  rotated  180 
degrees  (tapered  leading  edge  and  elliptical  trailing  edge). 
No  difference  was  noticed  in  tlie  results  for  lift  and  drag. 
Laitone  (1996  and  1997)  did  notice  a significant  increase 
in  lift  at  R.e„  = 20, 700  for  a thicker  reversed  NACA 
0012  airfoil  (the  sharp  trailing  edge  was  facing  the  flow). 
Further  tests  will  be  conducted  in  the  wind  tunnel  at  larger 
Reynolds  numbers  on  flat-plate  and  cambered-plate  wings 
to  complete  this  study  of  tlie  leading-edge  geometry  effect 


a (degrees) 

Figure  61:  Trailing-edge  geometiy  effect  on  drag 
coefficient  at  Re,,  — 80,000  on  flat  plates  in  the 
wind  tunnel  with  UND-FBl 
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Figure  62:  Trailing-edge  geometiy  effect  on  pitching 
moment  coefficient  at  Rcr^  ~ 80,000  on  flat  plates 
in  the  wind  tunnel  with  UND-FBl 
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Figure  63;  Trailing-edge  geometry  effect  on  lift 
coefficient  at  Rcc  ~ 80, 000  on  cambered  plates  in  the 
wind  tunnel  with  UND-FBl 
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Figure  64:  Trailing-edge  geometry  effect  on  drag 
coefficient  at  Re^  — 80, 000  on  cambered  plates  in  the 
wind  tunnel  with  UND-FBl 
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Figure  65:  Trailing-edge  geometry  effect  on  pitching 
moment  coefficient  at  Rcc  — 80, 000  on  cambered  plates 
in  the  wind  tunnel  with  UND-FBl 


on  the  aerodynamic  characteristics  of  thin  wings/airfoils. 

Effect  of  freestream  turbulence 

Mueller  et  al.  (1983)  showed  that  an  increase  in 
freestream  turbulence  intensity  reduced  the  minimum  drag 
acting  on  an  11%  thick  Lissaman  7769  airfoil  at  Rcc  ^ 
150,000  and  slightly  increased  This  was  caused 

by  an  earlier  laminar  shear  layer  transition,  hence  earlier 
flow  reattachment  (i.e.,  a shorter  separation  bubble),  with 
a larger  turbulence  intensity.  At  large  angles  of  attack 
where  the  flow  is  mostly  separated,  they  observed  an  in- 
crease in  drag  coefficient  with  an  increase  in  turbulence 
intensity.  Increasing  the  turbulence  intensity  also  helped 
to  eliminate  some  of  the  hysteresis  encountered  in  Ci  and 
Cd  for  that  particular  airfoil. 

Pohlen  (1983)  also  looked  at  the  influence  of  turbu- 
lence intensity  on  a 13%  thick  Miley  airfoil  (M06-13-128) 
(see  Miley,  1972).  He  found  that  increasing  the  turbu- 
lence intensity  helped  to  reduce  the  hysteresis  in  Q and 
Cd  and  slightly  improved  airfoil  performance. 

Tests  were  then  conducted  in  the  wind  tunnel  with 
different  screens  upstream  of  the  flat-plate  sAR  =^1.5 
models  and  a flow  restrictor  downstream  of  the  model  to 
see  if  a difference  in  the  turbulence  intensity  could  result 
in  different  aerodynamic  properties  for  the  models  used 
in  this  investigation.  The  flow  restrictor,  or  strawbox, 
was  made  of  drinking  straws  packed  in  a wooden  frame 
and  placed  between  the  test  section  and  the  diffuser.  The 
additional  turbulence  intensity  generated  by  the  strawbox 
was  determined  to  be  approximately  0.05%  (Brendel  and 
Huber,  1984).  Table  3 indicates  the  mesh  size  and  nomi- 
nal freestream  turbulence  intensity  in  the  test  section  with 
only  a screen  present  (no  flow  restrictor). 


Screen 

Mesh  size 
(meshes/ cm) 

Wire 

diameter 

(mm) 

Turbulence 

% 

Fine 

7.09 

0.245 

0.25 

Medium 

3.15 

0.508 

0.45 

Coarse 

0.64 

1.397 

1.3 

Table  3:  Turbulence  screen  data 
(Pohlen,  1983;  and  Brendel  and  Huber,  1984) 


No  measurable  differences  were  observed  in  the  re- 
sults for  different  turbulence  intensities  at  Re^  = 60, 000 
on  the  sAR  — 1.5  flat-plate  model,  as  shown  in  Fig- 
ures 66  and  67.  Only  a slight  increase  in  CL„,,^,,and  an 
increase  in  Cd  for  large  angles  of  attack  was  noticed 
for  the  case  with  the  fine  mesh  and  with  the  strawbox. 
All  other  cases  gave  the  same  results.  Therefore,  the  ef- 
fect of  turbulence  intensity  appeared  to  be  minimal  in  the 
wind  tunnel  for  the  models  tested  Similar  results  were 
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Figure  66:  Freestream  turbulence  effect  on  lift  coefficient 
at  Rcc  = 60, 000  in  the  wind  tunnel  for  the  sAR  =1.5 
flat-plate  model  with  UND-FBl 


Figure  68:  Schematic  of  corner  flow  on  wing 
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Figure  67:  Freestream  turbulence  effect  on  drag 
coefficient  at  Re^  ~ 60, 000  in  the  wind  tunnel  for  the 
sAR  = I Xj  flat-plate  model  with  UND-FBl 


obtained  at  Rxtc  ~ 120,000  in  tlie  wind  tunnel,  and  at 
Rcc  — 39,  000  and  Re^  = 60, 000  in  the  water  tumiel. 

Effect  of  endplates  on  2D  measurements 

It  has  been  shown  in  previous  experiments  at  Notre 
Dame  that  the  presence  of  the  endplates  during  2D  tests 
usually  leads  to  a larger  For  an  18%  thick  airfoil 

(NACA  663  — 018),  Mueller  and  Jansen  (1982)  showed 
that  the  interaction  between  the  endplates  and  die  model 
resulted  in  a 20%  increase  in  at  Reynolds  numbas 
between  60, 000  and  200, 000. 

Since  most  of  the  tests  in  the  current  investiga- 
tion were  conducted  at  very  low  speeds,  die  interaction 
between  the  boundary  layer  growing  on  the  endplates 
and  the  wing  created  a corner  flow,  as  depicted  in  Fig- 
ure 68,  which  acted  over  a significant  portion  of  the  wing 
span  and  significantly  altered  the  2-dimensionality  of  the 


flow  over  the  wing.  This  phenommon  of  die  comer 
flow  has  been  investigated  by  several  authors,  including 
Hawthorne  (1954)  and  Barber  (1978)  who  looked  at  the 
flow  around  stmts  near  a wall. 

In  order  to  verify  the  effect  of  the  endplates  on  the 
aerodynamic  characteristics  of  die  Eppler  61  airfoil,  a 
3-piece  Eppler  61  model  was  used.  Widi  tliis  setup,  a 
section  of  an  Eppler  61  model  was  free  to  move  between 
two  odier  sections  of  the  same  airfoil.  These  two  other 
sections  were  fixed  to  die  endplates  in  the  wind  tunnel, 
as  shown  in  Figure  69,  at  die  same  angle  of  attack  as 
the  middle  section.  A small  gap  was  present  between  the 
end  models  and  die  center  piece  connected  to  the  force 
balance.  Figure  70  shows  the  tliree  pieces  of  the  3-piece 
Eppler  61  model. 

The  angle  of  attack  of  the  3-piece  Eppler  61  model 
was  adjusted  to  a certain  value  and  the  velocity  in  the 
tunnel,  hence  die  Reynolds  number,  was  varied.  The  be- 
havior of  Cd  and  Q was  measured.  From  the  previous 
2D  results  on  die  Eppler  61  airfoil,  it  was  detemiined 
that  occurred  at  a = 0*^  and  the  angle  for  zero  lift 
was  about  a = —2'^.  The  behavior  of  Cd  vs  Rec  was 
first  obtained  at  diese  two  angles  of  attack.  Figures  71 
and  72  show  that  die  drag  coefficient  widi  die  3-piece 
Eppler  61  model  was  much  smaller  dian  with  the  full 
model.  This  result  is  similar  to  diat  reported  by  Mueller 
and  Jansen  (1982)  for  the  NACA  663  — 018  airfoil.  The 
lift  coeflicient  widi  the  3-piece  model  was  higlier  than 
with  the  full  model.  The  aerodynamic  characteristics  with 
the  3 -piece  model  were  closer  to  tme  2-dimensional  re- 
sults, where  a larger  Q and  smaller  Cd  would  nonnally 
be  expected.  The  behavior  of  Ci  and  Cd  with  Reynolds 
numbers  also  followed  die  expected  trends.  A reduction 
in  Cd  and  an  increase  in  Ci  were  observed  with  increas- 
ing Reynolds  numbas.  Results  from  Aldiaus  (1980)  and 
de  Vries  et  al.  (1980)  are  also  included  in  the  figures 
for  comparison.  As  was  mentioned  earlier  for  Altliaus, 
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Figure  69:  New-balance  arrangement  for  3-piece 
Eppler  61  airfoil  tests  in  wind  tunnel 


these  investigators  used  a strain  gauge  force  balance  to 
measure  lift  and  a wake  rake  to  measure  drag.  Since  the 
drag  measured  with  a wake  rake  is  usually  obtained  at 
the  mid-span  of  the  model,  it  does  not  take  end  effects, 
or  3D  effects,  into  account.  These  end  effects  can  be  sig- 
nificant at  very  low  Reynolds  numbers.  Therefore,  drag 
coefficient  results  from  Althaus  and  de  Vries  et  al.  were 
expected  to  be  smaller  than  the  present  results  and  this 
trend  was  obsen^ed. 

In  order  to  study  the  effect  of  endplates  at  low 
Reynolds  numbers,  Selig  et  al.  (1995)  showed  how  Cd  can 
var\^  for  a two-dimensional  airfoil  along  the  span  of  the 
model  at  low  Reynolds  numbers.  Figures  73  through  76 
show  drag  polars  on  an  SD6060  airfoil  at  Rcc  = 60, 000, 
100,000,  200,000  and  300,000,  as  obtained  using  a 
wake  rake  with  the  momentum  technique  for  Cd  and  a 
strain  gauge  force  balance  for  Ci.  For  Rcc  = 60,000 
and  Rcc  = 100, 000,  the  drag  polars  varied  significantly 
along  the  span  of  the  model,  which  implied  a three- 
dimensional  flow.  At  Rcc  = 200, 000  and  especially  at 
Rcc  = 300, 000,  the  drag  polars  were  relatively  constant 
along  the  span  and  a nearly  two-dimensional  flow  was 
believed  to  exist. 


(a)  Top  view 


Lower  piece  , 


Sting  hole 


(b)  Individual  pieces 


Figure  70:  3-piece  Eppler  61  airfoil  model  tested 
in  wind  tunnel 
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Figure  71:  Eudplates  effect  on  2D  characteristics  of 
Eppler  61  airfoil  at  a =■-  0°  wFh  UND-FB2 


Figure  72:  Endplates  effect  on  2D  characteristics  of 
Eppler  61  airfoil  a/  a ~ — 2°  with  UND-FB2 
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Figure  73:  Drag  polar  for  SD6060  airfoil  at 
Re^  ^ 60,  000  (Selig  et  al,  1995) 


Figure  75:  Drag  polar  for  SD6060  airfoil  at 
Rcc  = 200, 000  (Selig  et  al,  1995) 
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Figure  74:  Drag  polar  for  SD6060  airfoil  at 
Rcc  = 100, 000  (Selig  et  al,  1995) 


Figure  76:  Drag  polar  for  SD6060  airfoil  at 
Rcc  — 300, 000  (Selig  et  al,  1995) 
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(a)  Oil  flow  on  Clark  Y airfoil  at  Re,:  = 2.5  x 10'^ 
(Peake  and  Tobak.  1982) 


(b)  Conjectures  pattern  of  skin-friction  lines 
(Peake  and  Tobak.  1982) 


Figure  77:  Flow  pattern  over  a rectangular  wing  of 
aspect  ratio  3.5  at  angle  of  attack  beyond  stall 
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For  semi-infinite  models,  and  also  for  finite  models, 
the  flow  is  higlily  three-dimensional,  as  ^^  as  shown  b>^ 
Williams  (1996)  for  low  aspect-ratio  wings  at  liigli  angles 
of  attack  and  Winkelman  and  Barlow  (1980).  whose  work 
was  l1cvie\^ed  by  Peake  and  Tobak  (1982).  Figure  77 
shows  the  flow^  pattern  o\  er  a finite  rectangular  ing  at 
an  angle  of  attack  beyond  stall. 

Other  results  for  the  3-piece  Eppler  61  model  indi- 
cated that  the  zero-lift  angle  of  attack  was  close  to  —2.5“ 
instead  of  a — 2'\  so  tests  were  conducted  at  the  new 
zero-lift  angle  of  attack,  as  shown  in  Figure  78. 

For  the  full  model,  it  was  detcmiined  that  (L/D)n,ax 
occurred  at  a 8 - . The  3 -piece  model  was  then  tested 
at  an  angle  of  attack  of  8" . Unfortunately,  the  lift  force 
acting  on  the  middle  section  of  the  model  became  large 
enough  at  larger  Reynolds  numbers  to  deflect  the  lou  er  tip 
of  the  section,  which  resulted  in  a poor  aligmiient  with  the 
fixed  bottom  part  of  the  model.  Some  three-dimensional 
effects  were  then  belie^'ed  to  be  created.  liich  yielded  re- 
sults tliat  did  not  follow^  a tvN  O-dimensional  trend:  the  lift 
coefficient  started  to  decrease  wiih  increasing  Rey  nolds 


Figure  78:  Endplat.es  effect  on  2D  characteristics  of 
Eppler  61  airfoil  at  a:  = —2.5°  with  UND-FB2 


number.  It  is  known  that  a semi-infinite  model,  where 
the  bottom  endplate  has  been  removed,  will  have  a lower 
lift  coefficient  than  an  infinite  model,  or  pseudo-infinite 
model.  Figure  79  shows  the  results  at  a --  8°.  It  is  in- 
teresting to  note  the  significant  increase  in  Cf  just  before 
Rdc  = 60,000.  Tliis  larger  lift  coefficient  at  Reynolds 
numbers  greater  than  60, 000  was  mentioned  earlier  when 
looking  at  the  behavior  of  (f  with  angle  of  attack  for  dif- 
ferent Reynolds  numbers. 

Conclusions 

It  has  been  shown  that  sensitive  equipment  must 
be  used  to  perfonn  aerodynamic  measurements  on  small 
models  at  vety  low  Reynolds  numbers.  To  acliieve  tliis 
goal,  a new  platfonn  aerodynamic  balance  was  designed 
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Figure  79:  Endplates  effect  on  2D  characteristics  of 
Eppler  61  airfoil  at  a with  UND-FB2 


and  built  at  the  University  of  Notre  Dame.  The  balance 
shows  good  sensitivity,  repeatability  and  accuracy.  Tech- 
niques to  reduce  noise  and  increase  the  accuracy  of  the 
results  in  the  water  tunnel  below  Reynolds  numb  as  of 
40, 000  are  still  being  investigated. 

Moreover,  it  has  been  shown  that  cambered-plate 
wings  with  4%  camba  offer  better  aerodynamic  charac- 
teristics than  flat-plate  wings  for  a given  Reynolds  num- 
ber. Reducing  the  Reynolds  number  can  lead  to  poor 
performance  due  to  the  large  reduction  in 

for  both  flat-  and  cambered-plate  wings. 

max 

The  turbulence  intensity  in  the  tunnel  and  the 
trailing-edge  geometry  have  been  shown  to  have  a very 
small  effect  on  the  measurements  of  lift,  drag  and  pitch- 
ing moment  on  the  thin  models  tested  in  this  investigation. 
However,  tlie  presence  of  endplates  can  affect  the  results 
for  lift  and  drag  due  to  the  interaction  between  the  bound- 


ary layer  growing  on  the  endplate  and  the  flow  around  the 
wing.  The  presence  of  endplates  at  low  Reynolds  num- 
bers reduces  the  2D  lift  and  increases  the  2D  drag  that 
would  normally  be  obtained  with  a truly  infinite  model, 
according  to  tests  on  the  Eppler  61  airfoil.  Furthermore, 
for  a given  model,  adding  one  endplate  leads  to  an  in- 
crease in  lift  compared  to  the  case  without  an  endplate. 
Testing  finite  wing  models  without  endplates  appears  to 
be  desirable  at  low  Reynolds  numbers  because  it  elim- 
inates endplate  effects.  The  full  effect  of  the  endplates 
is  still  under  investigation,  especially  for  cambered-plate 
wings. 

Finally,  trends  of  versus  Rcc  should  be  inves- 
tigated by  fixing  the  wing  at  the  angle  of  attack  yield- 
ing C£)^^^.^and  then  varying  the  wind  velocity,  hence  the 
Reynolds  number,  in  the  tunnel.  This  seems  to  help  to 
eliminate  scatter  in  the  results.  Tests  using  this  technique 
on  various  flat  plates  have  indicated  the  expected  trends 
of  a reduction  in  C'^j^^.^^and  an  increase  in  with  an 
increase  in  Reynolds  numb  a. 
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